Cooperative excitation and many-body interactions in a cold Rydberg gas 
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The dipole blockade of Rydberg excitations is a hallmark of the strong interactions between 
atoms in these high- lying quantum states P, Qj. One of the consequences of the dipole blockade 
£h ' is the suppression of fluctuations in the counting statistics of Rydberg excitations, of which some 

^ ' evidence has been found in previous experiments. Here we present experimental results on the 

dynamics and the counting statistics of Rydberg excitations of ultra-cold Rubidium atoms both 
£NJ ■ on and off resonance, which exhibit sub- and super-Poissonian counting statistics, respectively. We 

£N) ' compare our results with numerical simulations using a novel theoretical model based on Dicke states 

of Rydberg atoms including dipole-dipole interactions, finding good agreement between experiment 
' ^ ' ■ and theory. 

: 

(3JT), PACS numbers: 32.80.Ee,42.50.Ct,03.67.Lx 

+-> ' 

Atoms excited to high-lying quantum states, so-called Rydberg atoms, arc highly polarizable and, therefore, can 
interact strongly with each other at large distances 0, Q ■ The study of ultracold gases has opened up new avenues of 
Q-i the investigation of strong atomic interactions. A wealth of potential applications is associated with the manipulation 
' of ultracold atoms excited to Rydberg states interacting through strong van der Waals (vdW) or long-range dipolar 
interactions, ranging from studies of strongly correlated quantum systems to quantum information Q. A key signature 
of interactions between Rydberg atoms is the suppression of fluctuations in the number of excitations due to the dipole 
■ . blockade. While evidence of the dipole blockade has been observed in several experiments 0-0] , prior studies were 
unable to demonstrate sub-Poissonian behaviour with a statistically significant confidence level. Also, those studies 
were carried out for resonant excitation only, leaving open the question of how Rydberg excitations in the strongly 
^ ■ interacting regime evolve for finite detuning from resonance. 
I— ~~ 'j The dipole blockade is a hallmark of the strong interactions between Rydberg atoms in an ultra-cold atomic gas. 
When an ensemble of atoms is irradiated by laser light resonant with a Rydberg excitation (either using a single laser 
or a multi-step excitation scheme), due to the interaction between Rydberg atoms, the excitation of a particular atom 
\Q , can be suppressed by a neighbouring one that is already in a Rydberg state. The radius of influence of an atom in 
this sense is called the blockade radius, which depends on the interaction strength and the linedwidth of the Rydberg 
excitation. As a result of the dipole blockade the excitation dynamics of the atoms are strongly correlated, leading to 
■ a cut-off in the excitation when all the available blockade volumes in the sample have been exhausted [3fJ . Another 
\Q \ way of describing the phenomenon is to view the atoms within a blockade volume as a 'supcratom', meaning that 
rather than individual atoms, collective states of several atoms are excited In the regime where the size of the 
sample is larger than the blockade radius, quantum correlations are expected to interconnect the whole sample. The 
description of this fully correlated regime is complex and the role of the atomic correlations on the atomic observables 
has recently received much attention 0]. In a simple mean field- type model for the interacting Rydberg atoms [28| 
each Rydberg atom experiences only an average interaction energy, which corresponds to neglecting all quantum 
correlations in the system. In order to reproduce the observed fluctuations and spatial correlations of the Rydberg 
| excitations, however, more sophisticated approaches are needed [Iol - ll4j |. 

While the full dynamics of Rydberg excitations in a cold gas of thousands of atoms is difficult to access experimen- 
tally, a number of experiments have observed signatures of the dipole blockade [15j . such as the suppression of the 
excitation of a single atom in the presence of a nearby Rydberg atom @, @] , the spatial correlations between Rydberg 
excitations [lfjj and indications of a suppression of fluctuations in the counting statistics near a Forster resonance [1[ ■ 
In this Letter we present a comprehensive experimental and theoretical study of the counting statistics of Rydberg 
excitations of ultra-cold rubidium atoms in a magneto-optical trap (MOT). 

As a full quantum mechanical treatment containing the states of all the atoms in the MOT is not feasible, the 
results of our experiments are analysed using an original theoretical model based on the well established Dickc model 
of quantum optics. The Dicke model, originally introduced for describing cooperative spontaneous emission or super- 
radiance [13 and subradiance [jjjj . is here modified by including the vdW interactions between the Dicke collective 
states (DCS). Our approach leads to a manageable size of the basis set for the simulations, which is only on the 
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order of the number of Rydbcrg excitations expected, that is, around 10 — 20 for our parameters, to be compared to 
our experimental numbers of ultracold atoms between tens and hundreds of thousands. The collective Dicke states 
contain the full statistical information about the collective Rydberg excitation, not only the average quantities, and 
therefore allow calculation of all the moments of the excitation statistics, in particular the measured variance of the 
Rydberg excitations. Furthermore, including a coupling between Dicke states of different symmetry leads to a par- 
tially non-Markovian dynamics, which results in better agreement with experimental findings. 

In our experiments, rubidium atoms in a MOT are excited to Rydbcrg states with principal quantum number n be- 
tween 50 and 80 using a two-photon scheme and detected, after field ionization, by a channeltron charge multiplier (l9j 
(for details see Supplemental Material). The detuning of the resulting coherent excitation from the ground state to a 
highly excited Rydberg state can be varied by changing the frequency of the second step laser at around 1015 nm(thc 
first step laser at 421 nm is detuned by around a GHz from an intermediate resonance in order to avoid population of 
that state). In the limit of vanishing interactions between the Rydberg atoms each excitation to a Rydberg state is 
independent of all the others, leading to a Poissonian excitation process. When the interactions are strong, however, 
the excitation processes are highly correlated and one expects sub-Poissonian counting statistics. These two regimes 
can be quantified through the Mandel Q-parameter [2(| defined as 

where ((AN) ) and (N) are the variance and mean, respectively, for the number N of Rydberg excitations. Exper- 
imentally, we determine Q by repeating an excitation-detection cycle for a fixed set of parameters between 50 and 
100 times and extract the mean and the variance from that counting distribution. Clearly, for Poissonian processes 
Q = 0, while sub-Poissonian processes are characterized by a negative Mandel parameter Q < (as, for example, 
in the case of the photon arrival statistics of squeezed light), with Q = — 1 in the limit of complete suppression of 
fluctuations. When Q > one speaks of a super-Poissonian counting statistics, which can be the result of technical 
noise or of enhanced fluctuations intrinsic to the dynamics of the system (see discussion below) . For detection effi- 
ciencies r] < 100% , the observed value of the Mandel parameter is Qd = r\Q [1, @. This detected value is used in all 
experimental figures in this letter. The Supplemental Material reports our experimental verification that a Poissonian 
distribution Qd = is measured in the direct photoionization process, in which case no suppression of fluctuations is 
expected. 

In Fig. [U an example of the excitation dynamics to the 71D 5 / 2 state for atoms on resonance is shown. Initially, 
a rapid increase in the average detected ion number is concomitant with a decrease in Qd- After 0.5 /its the growth 
in the detected ion number slows down while Qd fluctuates around —0.4. Taking into account our detector efficiency 
of ss 40% [III, the measurements indicate that on resonance a value close to Q = — 1 is obtained. In the long time 
limit the effect of a varying Rabi frequency across the sample, due to both the distribution of atoms and the spa- 
tial laser profile and line- width, cause a dephasing of the excitations, leading to a steady growth in the number of 
Rydberg atoms excited. In order to obtain reliable statistics with a resonable number of experimental realizations, 
we performed numerical simulations to verify how quickly the mean and standard deviation converge for Poissonian 
and sub-Poissonian distributions. For example, to measure the Mandel parameter of a sub-Poissonian process with 
Q » —0.5 to within ±0.1 requires around 50 realizations. 

The resonant excitation dynamics can be easily visualized as a saturation of the maximum number of superatoms 
that fit into the volume of the MOT, leading to a limiting value for Q around —0.7 in a 'close-packing' geometry 
of blockade spheres. Away from resonance, however, the picture is more complicated. Figure [2] shows the mean of 
the detected ion number and of Qd as a function of the detuning from resonance for a fixed time. As expected, 
for zero detuning Qd is highly negative, i.e. close to —1. Away from resonance, however, the counting statistics 
quickly becomes super-Poissonian with Qd as large as 2 — 3. One reason for this may be the effect of laser frequency 
fluctuations on the time scale of the data acquisition periods (around 30 seconds). On either slope a small shift in 
frequency corresponds to a large difference in the number of detected ions, leading to a large variance. Our theoretical 
model discussed below, however, indicates that the measured positive value of the Mandel parameter off resonance is 
also intrinsic to the excitation dynamics of the system. 

For a better understanding of the physics underlying our experimental observations we need a theoretical treatment 
of our system that gives access to the full counting statistics, and hence to the Q-parameter, both on and off resonance. 
Although recent theoretical work has addressed the expected counting statistics for on-resonant excitation of collective 
Rydberg states [13, El: so far there have been no predictions for the off-resonant case. We describe our system as an 
ensemble of iVo indistinguishable atoms, for which the laser excitation results in a single wavefunction superposition of 
a two level system composed of the ground state \g) and the Rydberg state \R), driven by a laser with Rabi frequency 
fl and detuning S. Our approach is based on the DCS (l7ll2ll|. characterized by the cooperative number r defining the 
symmetry and the number TV of Rydberg excited atoms, |r, N). As the laser excitation preserves the DCS symmetry, 
only states with the same symmetry number r are excited when starting from an ultra-cold atomic sample described 
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FIG. 1: Rydberg excitation dynamics in an ultra-cold gas. The mean number of detected ions (No) and Qd are shown, in a 
and b, respectively, as a function of pulse duration for the resonant excitation of the 71D 5 / 2 state. The blue laser detuning from 
the intermediate 6P3/2 F=3 state is 1 GHz. The number of ground state atoms is 8 x 10 3 , the average density 1.2 x 10 10 cm" 3 
and the two-photon Rabi frequency is 40 kHz. The error bars are derived from 50 experimental realizations. The continuous 
lines are the prediction of our complete model, while the dashed line in a is the theoretical prediction of a model based on the 
Dicke symmetric states only. All quantities are rescaled by the detection efficiency. 



by a totally symmetric DCS. When vdW interactions are taken into account, states with the same symmetry and 
those with different symmetries can be coupled. 

Assuming for the initial atomic wave-function the fully symmetrical non-dcgcncratc DCS \r = Nq/2,N) = \N) s 
and neglecting the vdW coupling between symmetrical and antisymmetrical DCS, we write the atomic wavefunction 
as jM* (t)} = J2n°=o Qn (0 \N) s ■ The time evolution of the amplitudes is given by the following coupled equations 
(with H = 1): 



a 7v = iNSciN 

-i§ \y/(N -N) (N+l)a N+1 + ^N(N -N+l)a N ^ 



(2) 
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FIG. 2: Detuning-dependence of collective Rydberg excitations. The mean number of detected ions {Nd) and Qd are shown for 
a range of frequencies around the two-photon excitation of state 71D!/ 2 in a and b, respectively. The (Nd) data are well fitted 
by a top truncated Lorentzian lineshape, except in a frequency range around the maximum. The continuous and dashed lines 
are the predictions of our complete model, with the calculated average coupling strength and 100 times that coupling strength, 
respectively. The black dot-dashed line simulates the effect of the laser jitter (assumed to be 0.5 MHz) on Qd, assuming that 
Q = independently of the detuning, based on the data of Fig. 2a. The blue laser detuning from the intermediate 6P 3 / 2 F=3 
state is 1 GHz. The number of ground state atoms is 5 x 10 4 , the average density 4 x 10 10 cm -3 , the two-photon Rabi frequency 
40 kHz and the excitation pulse duration 0.3 /is. 



The WdV interaction W^ s = W N (N - 1) /(N (N - 1)) introduces an additional term 

- iW s N s a N (3) 

to the equation of motion for a at, where Wq is the sum over all the interaction energies between pairs of Rydberg 
atoms [23]. The number of excitations in the Dicke ladder is equal to the number of Rydberg excitations, that is, 
a few tens in our experiment. The solution of the above equations allows one to calculate the cooperative mean and 
variance values, and all higher moments of the number of Rydberg excitations. It leads to the Rydberg excitation 
depicted by the dashed line of Fig. [TJ in reasonable agreement with the experimental observations for short times, 
but not in the long time limit. In addition it predicts asymmetric profiles for the laser detuning dependence in the 
blockade regime, not matching previous observations [23|] an d our data on the detected ions vs the laser detuning 
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shown in Fig. [2ji. 

The problem is that the vdW interaction couples the fully symmetrical DCSs to the remaining DCSs within a fixed 
iV subspace, which has a large degeneracy given by the binomial coefficient C^ Q . The most appropriate orthonormal 
basis within that subspace is composed of the symmetrical DCSs \N) and the superposition of non-symmetrical states 
{\N,q) ns }, with 1 < q < Cjy — 1, for which the vdW interaction matrix is diagonal. This original basis construction, 
detailed in the Supplemental Material, is far more general than the treatment of the vdW interactions provided we 
know the basis diagonalizing the coupling. Within this basis the vdW matrix contains the terms Wj] , WE, W^ q , see 
Supplemental Material. 

By writing the enlarged atomic wave-function as a superposition of symmetric and non-symmetric states, we obtain 
that the time evolution of Eq. is completed by terms describing the vdW coupling to the amplitude b q N of 
the non-symmetrical state. As a first approximation, we neglect the laser excitation of the amplitudes b q N because 
of their large degeneracy and the resulting weak occupation for each of them. After the elimination of the equations 
for bjy, the temporal evolution of the symmetric states determined by the coupling with a bath, made up of the the 
remaining collective states, therefore contains both Markovian and non-Markovian memory function contributions. 
For our experimental conditions, this quite general result can be simplified when calculating the memory functions see 
Supplementary Material. In the end for each number N of Rydberg excitations the evolution is described by a system 
of two equations, one for and one for a single level c/v defined as the superposition of all asymmetric DCS's of 
the fixed N subspace. Thus the bath is modeled by a number of states equal to that of the symmetric DCS's and the 
total required basis set required for the simulations is approximately twice the number of Rydberg excitations. The 
bath is Markovian and produces a shift of the symmetric DCS energy precisely compensating the vdW eigenenergy of 
Eq. ([3]). This compensation leads to symmetric excitation profiles as in Fig. [2] The vdW interaction for the symmetric 
DCS finally replaces Eq. d3]) by the following one: 

a N = -iW s N s c N , (4) 

For the non-symmetrical DCS 

c N = iSNc N - iW^ s aN- (5) 

Notice that within the above equations the vdW interaction acts as a coherent coupling between symmetrical and 
antisymmetrical states, leading to a variety of coherent evolution features, for instance the revivals appearing as nar- 
row peaks in the Q dependence on parameters as the laser detuning, see Fig. [2] or the laser pulse duration. However 
for each TV manifold the revival conditions are slightly different, and for a large number of Rydberg excitations, the 
revivals are washed out. The last step of this theoretical analysis is the introduction of the laser excitation into 
evolution of c/y. We treat it perturbatively, with a limitation on the temporal evolution imposed by the spectral 
linewidth Auj of the excitation laser. The final equations are reported within the Supplemental Material. 

Typical results of numerical simulations using the Dickc model are shown in Figs. Q] and [2] together with our 
experimental data. Clearly, the salient features of the experiment are reproduced well by our model, both in the dy- 
namics and in the dependence on the detuning. In particular, the off-resonant behaviour of the Q-factor, which in the 
experimental data features two symmetric peaks, is also visible in the simulation. The variation in the interatomic 
interaction curves, see Supplemental Material, is reflected in the uncertainty in our determination of the coupling 
strength for the theoretical curves in Fig. [2] Perfect quantitative agreement between experiment and theory is not 
expected as the Dicke model does not, at this time, include the spatially varying atomic density of the MOT, Dopplcr 
broadening, the spatial inhomogeneity of the excitation lasers, the full manifold of sublevels of the Rydberg states, 
and the effects of magnetic and possible residual electric fields, all of which might contribute to the larger linewidth 
observed experimentally. In particular, the dip in the Rydberg number around zero detuning seen in Fig[2] (a) could 
be smoothed out when taking into account the density variation over the cloud. In fact, Figure 3 shows that the 
variation of Q with the detuning depends strongly on the density (or, equivalently, on the C§ coefficient that describes 
the van-der-Waals interaction between Rydberg atoms). As the density is increased, the value of Q on resonance goes 
from zero to a limiting value of around —0.7, whereas the position and height of the off-resonant peaks symmetric 
around zero detuning increase steadily. In order to get better agreement, one could either perform experiments with 
close to uniform atomic densities (e.g., using flat-top excitation beams much smaller than the atomic sample) or 
develop a theoretical description that allows for the non-uniform density to be included. 

In summary, we have shown clear evidence of sub-Poissonian counting statistics of Rydberg excitations in an ultra- 
cold atomic sample, which are a clear signature of many-body quantum correlations in such a system. We have 
characterized both the dynamics and the dependence on the detuning of the Q-P aram cter and reproduced the main 
features of both using a novel model. In future studies, it will be interesting to look directly at the spatial correlations 
between Rydberg excitations at the counting statistics for the antiblockade regime of (24[, and at other, more 
indirect, signatures of Rydberg- Rydberg interactions, for example in EIT experiments (25j . 
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FIG. 3: Predictions for the detuning-dependence of collective Rydberg excitations, (a) Theoretical dependence of the Q- 
parameter on the laser detuning for a range of atomic densities (in cm" 3 ): 1 x 10 10 (blue), 4 x 10 10 (green), 10 x 10 10 (orange) 
and 40 x 10 10 (red), respectively . (b) Dependence on the density of the on-resonant Q (blue line) and the value of Q (red line) 
and position (black dashed line, right-hand scale) of the first off-resonant peaks. 
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Supplemental Material 

Sec. I provides additional information on the experimental setup, on the reliability of the measurements yielding 
negative values for Q d and on the absence of ions within the laser excitation volume. Within the theoretical description 
based on the DCS, after a brief recall of the mean-field Sec. II describes the construction of a basis of non-symmetric 
states the construction of the bath of the non-symmetric states and their excitation by the laser. Sec. Ill discusses 
the vdW parameters determining the Rydberg interactions. 



I. SET-UP 



Experiments are performed using a two-photon excitation scheme for Rydberg excitations in a small cloud of around 
10 5 87 Rb atoms trapped in a MOT of radius around 30 /xm and densities between 1 x 10 10 and 5 x 10 10 atoms/cm 3 
determined to within a factor of 2 from fluorescence images (for details see [Hj]). The sizes of the excitation beams 
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(around 200 /im) are chosen such as to obtain intense beams that arc approximately uniform across the atomic sample. 
Single atom two-photon Rabi frequencies of around 100 kHz are achieved. The excitation lasers are pulsed for up to 
20 /is, after which the Rydberg atoms are field ionized and the ions are detected by a channeltron. For the two-photon 
excitation, a laser beam at 421 nm is detuned around 0.5-1 GHz from the 5Si/2(F = 2) — > 6P 3/ / 2 (F' = 3) transition, and 
an infra-red beam of wavelength 1013-1015 nm provides the second step to the Rybderg state 6P 3 / 2 (F' = 3) — > nS/nD. 
Both of these beams are controlled using acousto-optic modulators (AOMs) with a rise-time of around 80 ns. Both 
lasers are stabilised to a scanning Fabry-Perot cavity to prevent long-term frequency drift. Wc estimate that the two- 
photon lincwidth on the timescale of the excitation pulses is around 300 kHz. In order to achieve reliable statistics, 
the excitation sequence is repeated up to a few hundred times with the repetition rate limited to 10 Hz in order to 
avoid a buildup of charge on the quartz cell surrounded by the external field plates [26j . The complete acquisition 
procedure takes around 30 s. The detection efficiency for ions in our experiment is around 40%. 

In order to verify that the observed sub-Poissonian statistics for resonant excitation is actually due to collective 
Rydberg excitations and not an artefact of the limitations of the experimental apparatus, we performed a series of 
additional experiments. Most importantly, we performed an experiment that interpolated between Poissonian and 
sub-Poissonian statistics for a fixed average number of Rydberg excitations, thus ruling out a spurious effect due to 
a saturation of the channeltron. To that end, the volume of the MOT was varied using a combination of magnetic 
field gradients, trapping beam sizes and laser powers chosen such as to keep the number of excited Rydberg atoms 
constant. As expected, for a small cloud the Rydberg dynamics was in the highly collective regime and a negative Qd 
was measured. As the volume increased with the number of Rydberg excitations kept fixed, the excitation dynamics 
tended to a Poissonian distribution with Q ~ 0. Moreover, we performed experiments in which the first step laser 
was tuned to resonance with the 5S!/ 2 (F = 2) — > 6P 3 / 2 (F' — 3) transition so that atoms were directly ionized through 
the absorption of two photons at 421 nm. Up until counting saturation, this process was found to obey Poissonian 
statistics with Qd = 0. 

The creation of spurious ions, that is, those not created by field ionization of Rydberg atoms, is a delicate issue 
in studies of the dipole blockade. The Stark shift due to a single ion is larger than the van-der-Waals interaction 
between two atoms in the investigated D 5 / 2 state. Thus the creation of one ion in the sample can lead to a blockade 
due to Coulomb interactions which suppresses the excitation of other Rydberg atoms and thus mimicks the dipole 
blockade. A possible cause for the presence of spurious ions is the absorption of two 421 nm photons leading to direct 
ionization from the intermediate 6 2 P 3 / 2 state [26|. By detuning the blue laser 1 GHz from this state excitation the 
probability of creating an ion is 3 x 10 -7 per 2 /is pulse. For the low atoms numbers in our experiment this equates to 
around 3 x 10 -2 per shot and thus is not a problem. Also, black-body radiation limits the lifetime of Rydberg atoms, 
for the 71-D5/ 2 state the effective lifetime is ~150 /xs (27j : however, the excitation and ionization is performed within 
< 10 fis and thus should have negligible impact on the counting statistics. Finally, at the typical atomic densities 
used in the experiment the decay process is not affected by super-radiance. Any loss from the initially prepared state 
will, therefore, follow Poissonian statistics. 

II. THEORY 
A. Mean field 

In a many-body system with interactions, the mean-field theory replaces all interatomic interactions with an average 
or effective interaction, reducing the many-body problem to an effective one-body problem. Within this approach the 
Rydberg state energy is shifted by the vdW interaction characterized by the parameter Wo given by a sum over all 
the Rydberg pairs inside the atomic cloud, see [28j . For instance, in a medium with uniform atomic density n and for 
a Cq interaction depending as r -6 on the interatomic distance r, Wq is approximatively J^^s ~ (-j-) n 2 . 

B. vdW basis in the cooperative approach 

The vdW interaction couples the fully symmetrical DCSs to the remaining DCSs within the fixed TV subspace 
having the large C^ o degeneracy. The most appropriate orthonormal basis within that subspace is composed of the 
symmetrical DCSs \N) S and the superposition of non-symmetrical states {\N, q) ns }, with 1 < q < C$ o — 1, for which 
the vdW interaction matrix is diagonal. To obtain the correct number of states within the C^ o ensemble, we first 
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remove from the {\N,q)} set one state labelled \N,p), assuming W^L = W^ s . Considering the state 

1 



\N,q) = \N,q) + 



°iV 



1 



\N,p), 



(6) 



for 1 < q < C$ — 1. Then the basis is constructed by subtracting for each |iV, q) state its projection on the symmetrical 
\N) state. We obtain the {\N, q) ns } states defined as 



\N,q) ns = [\N,q)- s (N\N,$)\N) s 



(7) 



It is easy to verify that such a basis {|iV) s , {|A^, <z) rls }} is orthonormal. The van der Waals coupling is diagonal with 
eigenvalue for the restricted basis {\N, q) ns }; these eigenvalues are not modified by the above basis changes. The 
coupling between the states \N) S and \N,q) na is 



W. 



N 



(N\W\N,q) 7 



w - w - 



Nl 



(8) 



The above construction of the basis is very general, provided we know the basis which diagonalizes the coupling. 
The only approximation is the existence of the state \N,p). For a large number of atoms, such an approximation is 
reasonable because it is possible to identify a state \N,p) such that W~ « W^ s . 



C. Enlarged wavefunction evolution 



By writing the enlarged atomic wave-function as 



\^(t)) = J2^N (t)\N) s +Y,b q N (t)\N,q) T 



(9) 



N=0 



the ajv time evolution of Eqs. (2) and (3) in the main text is completed by terms describing the (b N , a^) vdW coupling. 
As a first approximation, we neglect the laser excitation of the amplitudes b q N because of their large degeneracy and 
the resulting weak occupation for each of them. Thus eliminating the equations for b q N , we obtain the following 
integro-diffcrcntial equations: 

n r 



a N = -i(W* - NS)a N - i- \ y/(N -N) (N + l)a N+1 + (N Q - N + l)a JV _ 1 



+ E 



dr 2 dr ss 



W, 



i(W N -NS) 
o V gg ) 



-ajv (t — t) dr. 



The coupling &n — > appearing in Eq. (fit)]) is described through the correlation function /at(t) as 



djv = - / f N (r)e iNST a N {t-T)dT, 
Jo 

with the correlation function given by 



q °JV 



*L _ 2i4-W t N 
dr 2 dr - 



-iff T 



(10) 



(11) 



(12) 



where W^i was eliminated using Eq. (|5J|. With only two Rydberg excited atoms within the volume V, fN becomes 



Mr) 



wl 



2V/N 



N 



N 2 



v 2y2 iV 4 



■ 2w o 

e v dv. 



(13) 



The above integral has an analytical solution. The solution for the correlation function is composed of a long-memory 
component f l N (r) and a short-memory one, f N (r). The real and imaginary parts of f N (r) are represented in Fig. [4] 
The f N (t) contribution treated in the Markov approximation leads to the following imaginary term: 



- / fk ( T ) ex P (i5jT)a N (t - r) dr = iW^ s a N (t) 
Jo 



(14) 
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FIG. 4: Real and imaginary part of (t) vs the delay time r for typical parameters investigated in our experiment. 



This term exactly compensates exactly the blockade energy term introduced in Eq. (|10j) above for the symmetric DCS 
theory and, as a consequence, produces the symmetrical resonance excitation profiles reported in the Fig. 1 of the 
main text. 

The long-memory part can be treated exactly, but a good approximation, up to a few tens of Rydberg excitations, 

is 



JV\2 



(15) 



Using the above short-memory and long-memory parts, the vdW coupling can be described by representing the ensem- 
ble of the non-symmetrical levels with j Rydberg excitation in terms of a single level Cj defined as the superposition all 
asymmetric DCS's of the fixed N subspace. Within this simplified scheme for each number N of Rydberg excitations 
the evolution is described by the following system of two equations, one for and one for a single level cjv, reported 
also in the text: 



a N 



-iNSaisr — i 



y/(N -N)(N + l)a j+1 + (N - N + l)a N _ 



iW s N s c N , 



-iNScn - iW^ s a N 



(16) 
(17) 



D. Laser excitation of the non-symmetric bath 



The last step of this theoretical analysis is the introduction of the laser excitation into evolution of cat. We treat 
it perturbatively, with a limitation on the temporal evolution imposed by the spectral lincwidth Aw of the excitation 
laser. The final result is a modification of Eq. (|T7)) into 



cn 



— +iN5 I cjv - iW ss a N + — \c N -i\ , 



with all coherence terms cnc* n _ 1 set to zero. Tn depends on the laser linewidth Alo as 



n 2 (N - 2N) 



Aw/2 



Aw/2 



(Aw/2) 2 + (5 - NW SS ) 2 (Aw/2) 2 + {5 + NW SS ) 2 



(18) 



(19) 



III. RYDBERG INTERACTION PARAMETERS 



The numerical results rely on the precise values for the vdW coupling strengths. In a first approximation, we can 
consider that, in zero field, two Rydberg atoms experience a Ce/r 6 interaction, Cq being calculated on the basis of 
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perturbation theory of [31J . This model fails to describe the interaction accurately, first at small interatomic distances 
and secondly because of the Zeeman degeneracy of the two-atom Rydberg state. While at large interatomic distances, 
the two Rydberg states are dipole coupled via one additional state, at very short interatomic distances a large number 
of states is coupled, and the overall coupling is reduced. In our calculation of the interaction energy Wo, we introduce 
an appropriate cut-off radius (roughly twice the atomic size) because the excitation of two Rydberg atoms having 
a distance smaller than the cut-off radius is very unlikely. At intermediate interatomic distances the perturbative 
treatment of [3l| is not fully valid because of crossing in the molecular levels. Then a full diagonalization of the 
Hamiltonian leads to a more complex dependence of the interatomic coupling on the interatomic distance, Ce/r 6 
at large distance and C3/V 3 , i.e. a dipole-dipole interaction, at smaller ones. Finally, Walker and Saffman [32| have 
shown the Rydberg interaction is strongly dependent of the relative orientation of magnetic momentum of the colliding 
atoms of the pair. For the of Ds states this leads to a variety of 21 possible interaction curves whose strengths vary 
by more than two orders of magnitude. The average presented in [32j allowed us to obtain an effective dipole-dipole 
interaction, i.e., a single interaction curve for the calculation of the parameter Wq. The variation in the interaction 
curves is reflected in the uncertainty in our determination of the coupling strength for the theoretical curves in Fig. 
1 of the main text . 
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